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Abstract 

A theoretical framework is presented to treat hadronic observables within 
analytic perturbative QCD beyond the leading order of the coupling 
and for more than one single large momentum scale. The approach 
generalizes and extends the pioneering work of Shirkov and Solovtsov 
on an analytic strong running coupling. Some applications to hadronic 
observables at the partonic level are also discussed. 



1 Homage to Igor Solovtsov 

I (NGS) met Igor for the first time many years ago in the beginning of the 
nineties during a physics conference in Dubna. It was a chance encounter 

1 Invited contribution in memory of I.L. Solovtsov at Seminar in Bogoliubov Laboratory 
of Theoretical Physics, JINR, 141980 Dubna, Russia, Jan. 17, 2008. 
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during the conference dinner, when it happened to get a seat next to him. 
Igor introduced himself and I took the opportunity to discuss with him about 
variational perturbation theory and its convergence properties, a subject on 
what he was working at that time, known to me from his publications. We 
soon found out that we had a lot of common interests in physics and we 
kept discussing for hours, while emptying a bottle of vodka. During this first 
encounter, I could, of course, not imagine that years later I would engross 
myself so strongly in Analytic Perturbation Theory, a subject pioneered by 
Igor and Dmitry V. Shirkov, as it actually happened. This activity will 
be surveyed in this invited contribution to his memory. Since then, I met 
Igor on several occasions, in Dubna and abroad, and I still have a strong 
recollection of our discussions. Igor indelibly forged his name in the annals 
of physics by his scientific achievements — no doubt. But those of us, who 
had the privilege to know him in person, will sadly remember and miss his 
intellectual creativity along with his kindness. 

2 Introduction 

Traditionally, perturbation theory in QCD suffers from an artificial singular- 
ity at momenta close to Aqcd, called (at the one- loop level) the Landau pole. 
Since the early days of QCD, theorists have tried different remedies to avert 
this problem, like infrared (IR) cutoffs, "freezing", etc. But approximately 
ten years ago, Shirkov and Solovtsov have devised an approach that avoids 
this problem by appealing only to a few basic principles of Quantum Field 
Theory — chief among them, Causality and Renormalizability — and avoiding 
the introduction of extraneous infrared (IR) regulators. Since this pioneering 
work appeared in the year 1996 [1, 2], the analytic approach to QCD pertur- 
bation theory has evolved and considerably progressed, shifting the cutting 
edge significantly (see [3-5] for reviews and further references), and finally 
culminating into the so-called Analytic Perturbation Theory (APT). Mean- 
while this analytic approach has been extended beyond the one-loop level 
[6, 7] and important techniques for numerical calculations have been devel- 
oped [8-13]. Also applications to the ultra-low momentum region have been 
carried out, e.g., [14], and alternative formulations of the strong coupling 
below the Landau pole have been proposed aiming to incorporate nonper- 
turbative input [15, 16]. 

The simple analytization concept of the strong running coupling has been 
generalized to the level of hadronic amplitudes [17, 18] and new techniques 
have been developed to deal with more than one large hard scale in the pro- 
cess [19,20], including also Sudakov resummation in exclusive processes. In 
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the later course of these investigations, it was realized that logarithms of 
the aforementioned second large scale — which can be the factorization or the 
evolution scale — correspond to non-integer (fractional) powers of the cou- 
pling, giving rise to Fractional Analytic Perturbation Theory (or FAPT for 
short) [21, 22]. At the heart of this development was the Karanikas-Stefanis 
(KS) [17] analytization principle which demands that all terms in a QCD 
amplitude that can affect the discontinuity across the cut along the negative 
real axis — oo < Q 2 < 0, and hence contribute to the spectral density, have 
to be included into the analytization procedure, i.e., the dispersion relation. 
The KS procedure encompasses the Shirkov-Solovtsov analytization concept 
of integer powers of the strong coupling and paves the way to the analyti- 
zation of any real power both in the Euclidean [21] as well as and in the 
Minkowski region [23]. The crucial advantages of this scheme are: 

• A diminished sensitivity on the factorization scale of typical QCD hard 
processes, like the factorized part of the pion's electromagnetic form 
factor — verified in [22] to the next-to-leading order (NLO). 

• A quasi renormalization-scale setting and scheme independence of the 
same observable at the NLO level [24] (see for an abridged version [25]). 

• The inclusion of evolution effects due to the running of the strong 
coupling beyond the one-loop order [23]. 

• A faster convergence of the perturbative expansion in terms of analytic 
images of real powers of the strong coupling [21]. 

• The resummation to all orders of the tt 2 terms, induced by analytic 
continuation into the timelike domain. This has been exemplarily ver- 
ified in [23] for the scalar Higgs decay into a bb pair at the four-loop 
level. 

The most important consequence of the analytic approach is that it ties 
the abstract mathematical requirements of causality and renormalizability 
to something tangible, like the calculation of the Bjorken [26] and the Gross- 
Llewellyn Smith sum rule [27], or the inclusive decay of a r-lepton into 
hadrons [28-30], the pion's electromagnetic form factor [19,20,22,24], the 
Higgs boson decay into a bb pair [23], and many other processes. In the 
present exposition we will present the bedrock of this approach, focusing our 
attention to selected applications and results beyond the leading order (LO) 
of perturbative QCD — conventional and analytic. 
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3 Basic Structure of FAPT 



Let us briefly review FAPT, specifying our notation, explaining its main 
principles, and describing its methodology. 
The running strong coupling in QCD 

a s {Q 2 ) = ^-a s {L) with L = In § , (1) 

where A denotes the characteristic scale of QCD, Aqcd, satisfies the 
renormalization-group (RG) equation 

with known /^-function coefficients up to the displayed order. [Their ex- 
plicit expressions can be found, for instance, in [23]]. At the one-loop order, 
a s (L) develops a Landau pole, while the two-loop solution of Eq. (2) has a 
square-root singularity, with more complicated singularities for still higher 
orders. Shirkov and Solovtsov [2] have shown that the ghost-singularity 
problem can be solved only on account of renormalizability — in terms of the 
RG equation — and causality — expressed in the form of a dispersion relation. 
Then, one obtains in Euclidean space analytic images of the coupling at loop 
order I 

A%(Q 2 )= \a^(Q 2 )] (3) 

L w J an 

following from the Kallen-Lehmann spectral representation 

with the spectral density at one loop given by (see [21, 23] for higher loops) 
p v [a) = - Im cr(-o-) = — —j^- (5) 

71" [7T 2 + L 2 (<t)] 

As a result, one then finds in the Euclidean space at one loop 

4 1} (Q 2 ) = i-p^r ( 6 ) 

while its counterpart in Minkowski space reads 

2lS 1} (s) = - arccos ( ] (7) 
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with L s = In (s/A 2 ). This procedure is sufficient in considering QCD pro- 
cesses with one large scale, but fails for non-integer powers of a s and is unable 
to accommodate terms, like 

• [a s (L)] 7o/2/3 ° < ► RG at one loop 

• [a s (L)] n ln[a s (L)] < — ► RG at two loops 

• [a s (L)] n L m < — ► Factorization 

• exp [— a s (L)F(x)] < ► Sudakov resummation (symbolically) 

which typically appear in perturbative calculations beyond the leading order 
due to the reasons already explained. Such terms do not modify the ghost 
singularities but they do contribute to the spectral density and, hence, their 
analytic images are inevitably required for the dispersion relation. In actual 
fact, precisely such terms are tantamount to fractional (real) powers of the 
strong coupling [21,22]. 

The core feature of FAPT is, as mentioned in the Introduction, the KS 
analytization principle [17]. The use of this principle allows the inclusion 
into the dispersion relations of logarithmic terms, of the sort ln(Q 2 ///p), 
or products of such logarithms with powers of the running coupling. To 
appreciate its meaning and usefulness, we present and compare different 
analytization concepts in Fig. 1. 

In this figure, the linear operations Ae and Am define, respectively, the 
analytic running couplings in the Euclidean (spacelike) 

n 



= A$ with AiHQ 2 )^ r^da (8) 

Jo (J -r v 

and the Minkowski (timelike) region 

:2l« with 2#>( a ) = r^lda. (9) 

Js °~ 



The above analytization operations can be represented by the following 
two integral transformations from the timelike region to the spacelike region 
(see, e.g., [4]): 

£«)] = 4P with ^)(Q 2 ) = Q 2 r 7 ^Ld* (10) 

Jo cr + Q 2 



and for the inverse transformation 



R[A®] = 2l« with ^) { s)^^J~^^-da. (11) 
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(a) | APT 

PT | PT | 

[ai l) (Q 2 )}" In (Q 2 / IV 2 ) 

I 1 I 1 

I 5a® (s) ^ ^'(Q 2 ) ; ; s$>(s) ^ £l°(Q 2 )j 

I R=D~ 1 j | R^D- 1 j 
I I I I 



(6) FAPT 



(c) FAPT 



Fig. 1: Illustration of different analytization concepts, (a) APT (b) FAPT, 
only strong-coupling powers, (c) FAPT, products of strong-coupling powers 
and logarithms. In APT, the index n is restricted to integer values only; in 
FAPT v can assume any real value. Further explanations are given in [23]. 



These two integral transformations are connected to each other by the rela- 
tion 

DR = RD = 1 , (12) 

valid for the whole set of analytic images of the powers of the coupling in the 
Euclidean as well as in the Minkowski space, {„4 n ,2l n }, respectively, and at 
any desired loop order of the perturbative expansion. 

In the spacelike region, the analytic images of the coupling can be ex- 
pressed in terms of the reduced transcendental Lerch function F(z,v) to read 
[21] (L = ln(Q 2 /A 2 )) 2 

1 F(e~ L 1 - u) 
ML) = ±- ' rw > ■ (13) 

where the first term corresponds to the conventional term of perturbative 
QCD and the second one is entailed by the pole remover (cf. l/(e L — 1) at one 



2 Everywhere in this presentation Greek labels denote non-integer (real) powers or in- 
dices. 
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loop) . This function is an entire function in the index v and has the properties 
A (L) = 1, A- rn {L) = L m for m G N, and A m (±oo) = for m > 2, m G N, 
while for \L\ < 2tt, it reads A U (L) = - [1/T(u)] J2Zo C(l - "~ [(-£) r A*] • 
In the timelike region, these images are completely determined by elementary 
functions [21] (L s = ln(s/A 2 )): 



sin 



(i/ - 1) arccos 



■k{v — 



(14) 



from which, for example, we get Sto(-^s) = 1 an d 2Li(£ s ) = L s . The salient 
characteristics of FAPT in comparison with APT and the standard pertur- 
bative expansion in QCD are compiled in Table 1, while for further reading 
and graphic illustrations we refer the reader to [21-23]. 



Table 1: FAPT versus APT and standard QCD perturbation theory (SPT) 



Theory 


SPT 


APT 


FAPT 


Space 




{A m } meN 




Series expansion 


Zf m a m (L) 

m 


EfmA m (L) 
m 


EfmA m (L) 
m 


Inverse powers 


[a(L)]- m 




A- m (L) = L m 


Multiplication 


a^a u = a^ +u 






Index derivative 


a? \n k a 







4 FAPT Applications at the NLO and Beyond 

In this section, we concentrate on applications of the presented FAPT formal- 
ism to two QCD processes beyond the leading order of perturbative QCD. 

The first process to be considered is the factorizable part of the pion's 
electromagnetic form factor at NLO accuracy in Euclidean space. This pro- 
cess has been widely discussed in the literature using various techniques — 
see e.g., [24] for a recent comprehensive analysis and comparison with the 
available experimental data. At leading twist, one has the convolution 
(A(z) ® zB(z) = Jo dzA{z)B{z)) 

Fj act (Q 2 ) = <p K (x, /4) ® T£ LO (x, y, Q 2 ; fi 2 F , fi 2 R ) ® ^(y, 4) , (15) 
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where the twist-two pion distribution amplitude (DA) (using x = 1 — x) 

<p w (x, /i 2 ) = Gxx \l + a 2 (p 2 ) Cl /2 {2x - 1) + a 4 (/U 2 ) cf 2 (2x - 1) + . . .1 (16) 

contains all non-perturbative information on the pion quark structure in 
terms of the Gegenbauer coefficients a n , determined at some typical hadronic 
scale « 1 GeV 2 [31-33]. Note that the quantity Fj act (Q 2 ) depends, be- 
yond the LO, on two scales: the factorization scale and the renormaliza- 
tion scale /xr. 

To appreciate the differences among the various analytization schemes, 
consider the scaled hard-scattering amplitude entering Eq. (15) in the so- 
called "Naive Analytization" (Naive- An) scheme [19, 20] in comparison with 
the "Maximal Analytization" (MA) scheme [24] (with the renormalization 
scale set equal to //^ = XrQ 2 , Ar being a numerical parameter): 



Naive Analytization 

I Naive— An 

A {2 \\ K Q 2 



[Q 2 T H (x,|/,Q 2 ; / u|,A R g 2 )] Naivc _ An = A? (XrQ 2 )^ (x,y) 



2 



i 1} (^y;AR,^) (ir) 



An 

Maximal Analytization 
[Q 2 T R (x,y,Q 2 ;4,\RQ 2 )] Max _ An = A {2) (\ R Q 2 ) t£\x,y) 

4 2 Wh(i) ' 



4ir 



4 1} (^2/;Ar,^) ■ (18) 



In these equations, t^\x,y) and t^ (x,y;\R, stand, respectively, for 
the LO and NLO hard-scattering amplitudes, computed in [34]. 

The Naive Analytization just replaces the strong coupling and its powers 
by their corresponding analytic images. This procedure is, strictly speak- 
ing, incorrect [20] because [„4i(L)] n / [a™(L)]An> owing to their distinct 
spectral representations. This scheme, whatever its theoretical shortcom- 
ings, works phenomenologically rather well [19,20]. Its direct improvement 
in [24] adopts instead the Maximal Analytization, which associates to the 
powers of the running coupling their own dispersive images, trading this 
way the usual power series expansion for a non-power functional expansion, 
i.e., [a™(L)]Max-An = A n (L). The difference between the two analytization 
schemes becomes apparent by comparing in the equations above the NLO 
terms. This theoretical improvement entails a phenomenological improve- 
ment as well. From Fig. 2 we see that the crucial advantage of the FAPT 
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analysis is that the dependence of the prediction for F^ act (Q 2 ) on the pertur- 
bative scheme and scale setting is diminished already at NLO. Next we will 
show that applying the KS analytization procedure, the result will become 
insensitive also to the variation of the factorization scale. 




2 3 4 5 6 7 8 12345678 

Q 2 [GeV 2 ] Q 2 [GeV 2 ] 



Fig. 2: (Left) Results for Q 2 F^ ct (Q 2 ) vs. Q 2 with /4 = Q 2 , fi 2 = 5.76 GeV 2 
in SPT of QCD (dashed line), using Naive Analytization (dash-dotted line), 
and with Maximal Analytization (solid line). (Right) The same quantity (in 
Max- An) in comparison with experimental data (see [24]). The broken lines 
denote the region accessible to the asymptotic pion DA, while the shaded 
strip marks the region of predictions derived with the pion DAs from nonlocal 
QCD sum rules [31] (cf. Eq. (16)). 



Consistent with this requirement the analytization of the logarithmic 
term ln(Q 2 //z F ) = ln(X R Q 2 / A 2 ) — ln(AR/x F /A 2 ) has to be performed as well, 
so that after some manipulations, explained in [22], we obtain 



[Q 2 T u (x,y,Q 2 ;4,X R Q 2 )] 



2M An 
KS 



A?(X R Q 2 )t^(x,y) 



2W(0), 



+ 



■4 2) (ArQ 2 K(i) 



4tt 



A? (XkQ 2 ) I 



+ 



47T 



,(0) 



C F t£>(x,y) (6 + 21n(xy)) , (19) 



with the deviation from the second line in Eq. (18) being encoded in the 
term 

A^(Q 2 ) ee 4 2 ) (Q 2 ) - 4 2 ) (Q 2 ) In [Q 2 /A 2 ] 

where 



(20) 



4 2) (Q 2 



a 



(2) 



(Q 1 



2 m^ 



2\ 1 An 



4vT 



a { P (Q 2 ) 



An 



KS 



(21) 
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Performing the analytization [22], we find 



4 2) (Q 2 ) 



Air 



4 2) (Q 2 )+c^fc(Q 2 ) 



(22) 



where 



h{Q 2 ) = Yl ma-n-i) 



n>0 L 




(23) 



and C(z) is the Riemann zeta- function. One can show (see for details [22])) 
that calculating F^ act (Q 2 ) under the proviso of the KS analytization provides 
an expression which is extremely stable against variations of the factorization 
scale. Indeed, varying the factorization scale from 1 GeV 2 to 10 GeV 2 , the 
form factor changes by a mere 1.5 percent and reaches just the level of 
about 2.5 percent for a (hypothetical) factorization scale of 50 GeV 2 . The 
sensitivity on the factorization scale using the Maximal Analytization is also 
a mild one, but the corresponding variation is, in round terms, two times 
larger. As regards the scale behavior of the form factor, both analytization 
schemes yield almost coincident results. Hence, [Q 2 F^ act (Q 2 )~j Kg in Fig. 
2 cannot be differentiated from [Q 2 F^ act (Q 2 )] Max _ An - In concluding this 




analysis, using FAPT the dependence on all perturbative scheme and scale 
settings, including the factorization (evolution) scale, is diminished already 
at the NLO level. 

We turn now to the second application, this time in Minkowski space: 
the decay of a scalar Higgs boson to a bb pair at the four-loop level of the 
quantity i?g from which one can obtain the width r(H — > bb) [23]. In this 
case, we will encounter no ghost singularities — in contrast to the Euclidean 
space. However, the analytic continuation from the spacelike to the timelike 
region will entail so-called "kinematical" it 2 terms, whose contribution may 
become with increasing order of the perturbative expansion as important as 
the expansion coefficients. 

To get a handle on the Higgs-boson decay, we consider the correlator of 
two scalar currents jf = ^>b^b for bottom quarks with mass m?,, coupled to 
the scalar Higgs boson with mass Mr and where Q 2 = —q 2 : 



Then, Rs{s) = ImII(- s — ie)/(2it s) and one can express the width in terms 
of R s , i.e., 





(24) 



r(H -► bb) = 



G F 



M n m 2 b (M n )R s (s = M 2 ) . 



(25) 
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One, finally, obtains Rs via the analytic continuation of the Adler function D 
into Euclidean space by applying on it the linear operation Am (equivalently, 
the integral transformation R), according to the analytization machinery 
illustrated in Fig. 1 This means that one has to calculate the quantity [36] 



Rs(s) = Rs(Q 2 = s,v 2 = s) = 3m 2 (s) 



1 + J> n <( S ) 

n>l 



(26) 



where the expansion coefficients r n contain characteristic tt 2 terms origi- 
nating from the integral transformation R of the powers of the logarithms 
appearing in Ds- The latter is related to Rs(s,s) by means of a dispersion 
relation. Notice that these logarithms have two different sources: one is the 
running of a s in D$ , while the other is related to the evolution of the heavy- 
quark mass m 2 (Q 2 ). As a result, the coefficients r n in (26) are connected to 
the coefficients d n in D$ (calculable in Euclidean space) and to a combina- 
tion of the mass anomalous dimension ji and the /3-function coefficients bj, 
multiplied by tt 2 powers [35-38]. 

The running mass in the l-loop approximation, mm, can be cast in terms 
of the renormalization-group invariant quantity rhm to read 



™fi){Q 2 ) = rn 2 (l) [a s (Q 2 )r f {l) (a s (Q 2 )) , 
where the expansion of fm(x) at the three-loop order is given by 



(27) 



f(i)(a s ) 



1 + a s 



26 



7i 
bi 



To 
bo 



+ a 



2 b\ 



S 16 6§ 



To 
bo 



Ti 



+ 2 



To 
bo 



7i 
h 



bob2 


( T2 


7o V 


b 2 i 


\bi 


bo). 



+ 0(a 3 s ) . 



(28) 



We are now ready to consider the analytization of the Adler function 



D s (Q 2 ;/j, 2 ) = 3mt(Q 



n>l 



(29) 



Expanding the running mass in a power series, according to 



m 2 {Q 2 ) = m 2 {a s {Q 2 )) 



m>l 



(30) 
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and choosing fi 2 = Q 2 , we find 



[3^]-; Dl\Q 2 ) = (a® (OP)) 1 * + J> (a(0(Q^ V ' +M ' 

n>l 

CO 

+ E^ ) (4' ) (Q 2 )) (3D 



m>l 



with 



min[i,m— 1] 

AW = e«+ £ *«2U- (32) 
fe>i 

Note that we have purportedly separated the mass-evolution effects (col- 
lected in the third term of Eq. (31)) from the original series expansion of 
D (truncated at n = I), the latter being represented by the second term on 
the RHS of Eq. (31). In practice, for Q > 2 GeV, i.e., for a s < 0.4, the 
truncation at m = I + 4 of the summation (30) produces a truncation error 
much smaller than 1 percent. 

Finally, we obtain i2g [FAPT from the quantity D<p(Q 2 ) by applying the 
analytization operation Am: 

jf MFAPT = Am[£) ( 0] 

(33) 



o - 2 

= 3m (J) 



n>l 

where we have used the short-hand notation 



■m+u 

n>l m>l 



Ws(s)\ n = a^\s)^(^-J^Hs) (34) 

and bo = tt-Ca — |Tr./V/ with Ca = N c = 3, Tr = |. The above expression 

contains, by means of the coefficients A$ ( ejjP) and the couplings fl^_„ , all 
renormalization-group terms contributing to this order, while the resummed 
7r 2 terms are integral parts of the analytic couplings by construction. 

The results for the quantity Rs(M^), calculated within different ap- 
proaches in the MS scheme, versus the Higgs mass Mjj are illustrated in 
Fig. 3. The long-dashed curve in this figure shows the predictions obtained 
by Baikov, Chetyrkin, and Kiihn [35] employing standard perturbative QCD 
at the I = 4 loop level of expansion. The solid curve next to it represents 
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60 80 100 120 140 60 80 100 120 140 



Fig. 3: Illustration of the calculation of the perturbative series of the 
quantity Rs(M^) in different approaches within the MS scheme: Stan- 
dard perturbative QCD [35, 36] at the loop level I = 4 (dashed red line 
with Atv /= 5 = 231 MeV), BKM estimates, by taking into account the 
O((a s y o ^ 4 (a s ))-terms, [37]— (dotted green line with Ajv /= 5 = HI MeV), 
and MFAPT from Eq. (33) for Nf = 5 (solid blue line), displayed for two 
different loop orders: I = 2 (left panel, Aj\r /= 5 = 263 MeV) and I = 3 (right 
panel with Ajy /= 5 = 261 MeV). The value of Ajv /=5 MeV in all cases corre- 
sponds to 2lf } (s = m 2 z ;N f = 5) = 0.120. 

the outcome of the FAPT machinery (cf. (33)), including in the second sum 
all evolution effects up to m = I + 4 and fixing the active flavor number to 
Nj = 5. Bear in mind that the ir 2 terms, induced through the analytic con- 
tinuation, are contained in the expansion coefficients o^ +UQ . On the other 
hand, the contributions of the higher-loop renormalization-group dependent 
terms are accumulated in the coefficients A„ by means of the parameters 
7i and bj. It is obvious that for this observable the standard perturbative 
QCD approach and FAPT yield similar predictions, starting with the two- 
loop running. The reason for the slightly larger FAPT prediction lies in the 
fact that the coefficients a u contain the resummed contribution of an infinite 
series of 7r 2 -terms that renders them ultimately smaller than the correspond- 
ing powers of the standard coupling. [The interested reader is referred to [23] 
for further details.] Finally, the lower (green dotted) curve by about 8% in 
both panels of Fig. 3 gives the estimate of Broadhurst, Kataev, and Maxwell 
(BKM) [37], which relies upon the so-called "naive non-Abelianization" and 
an optimized power-series expansion that makes use of the "contour integra- 
tion technique". 3 [Some more technical remarks can be found in [23], where 
the common elements between this approach and FAPT are worked out.] 

3 We thank A. L. Kataev for useful remarks pertaining to this figure in [23]. 
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5 Conclusions 



The generalized KS analyticity requirement [17] has proven successful in de- 
scribing hadronic observables at the partonic level for a variety of reactions. 
Although this requirement has to be extra postulated, it is the one that pro- 
vides a natural extension of the analyticity demand on the running coupling, 
proposed by Shirkov and Solovtsov [2], giving us a much broader under- 
standing of analytization. We have shown that including into the dispersion 
relations the contributions stemming from all terms that affect the spectral 
density (even though these terms do not influence the nature of the ghost 
singularities of the standard power-series perturbative expansion), makes it 
possible to treat processes containing two large momentum scales. Such ad- 
ditional scales, like the factorization or the evolution scale, enter in the form 
of typical logarithms whose incorporation into the spectral density naturally 
amounts to non-integer (fractional) powers of the coupling. This analyti- 
zation formalism — Fractional Analytic Perturbational Theory, developed in 
[21-23] on the theoretical basis of [17] — works equally well in both the space- 
like region (Euclidean space) as well as the timelike region (Minkowski space) . 
In the first case, the obtained expressions for the hadronic observables are 
singularity-free and turn out to be insensitive to the renormalization scheme 
and scale adopted, while bearing little sensitivity to the factorization scale, 
as well. In the timelike regime — where ghost singularities are absent — a bet- 
ter stability is achieved in terms of expansion coefficients up to a high loop 
order that in situ resum all it 2 terms induced by analytic continuation. 

In this short exposition we have not been exhaustive. We note in passing 
that we have derived closed-form expressions for the analytic-coupling images 
at the one-loop level in the spacelike [21] and in the timelike region, and fur- 
ther approximate expressions at the two- loop level [23], the latter supported 
by exact numerical results [11]. What is perhaps more, our approach provides 
a handle on the computation of power corrections and their coefficients to 
different hadronic reactions. We have already obtained leading-order power 
corrections to the pion's electromagnetic form factor and to the cross section 
of the Drell-Yan process into a lepton pair [17, 18], which put the developed 
scheme into a larger theoretical and phenomenological context. 
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